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Producing rare isotopes through statistical multifragmentation is investigated using the Mekjian 
method for exact solutions of the canonical ensemble. Both the initial fragmentation and the the 
sequential decay are modeled in such a way as to avoid Monte Carlo and thus provide yields for ar- 
bitrarily scarce fragments. The importance of sequential decay, exact particle-number conservation 
and the sensitivities to parameters such as density and temperature are explored. Recent measure- 
ments of isotope ratios from the fragmentation of different Sn isotopes are interpreted within this 
Oji picture. 
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I. INTRODUCTION 

> ■ 

The study of rare isotopes is attracting increasing attention due to the recent development of radioactive beam 

facilities where isotopes are produced through nuclear collisions. Projectile and target nuclei might vary in size from 
a few dozen nuclcons all the way to Uranium. The mechanism for rare isotope production might entail the transfer 
of a few nucleons between the projectile and target, induced fission of the projectile, or at highest energy transfer, 
multifragmentation. This last mechanism, which assumes temperatures of a few MeV, is the focus of this study. 
Multifragmentation is an extremely complicated process that takes place at time a time scale of 100 fm/c. At such 

"~^5 excitations tunneling through barriers plays a less central role than in fission. In this short time the system is able to 
sample an enormous number of configurations which makes dynamical calculations computationally very intensive. 

,_!^ ' However, this very complexity lends justification to statistical calculations which have proven remarkably successful 
in predicting mass yields for excitation energies in the range of a few MeV per nucleon jj],[| ■ 

It is our aim to apply simple statistical calculations to the study of the production of rare isotopes, nuclei with 
neutron-to-proton ratios far from the valley of stability. Our investigation addresses a variety of questions: 
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•rH . 1. How do yields depend on the physical parameters of the thermalized system? Such parameters are size, overall 

/\ ' neutron excess, density and temperature. 

h : 

2. How sensitive are the yields with respect to various aspects of the modeling such as breakup density, level 
densities, and barrier penetration probabilities? 

3. Is there a qualitative and possibly quantitative explanation for the isospin amplification effects seen for light 
isotopes? Furthermore, do these effects carry over into the production of heavier or exotic isotopes? It has been 
experimentally observed that the ratio of mirror nuclei, e.g. 15 0/ 15 N or t/ 3 He, can be of order 10 even though 
the neutron-to-proton ratio of the colliding nuclei is approximately 1.2 p|. 

4. By dividing the ratio of isotope yields from the fragmentation of 124 Sn by the same ratio from the fragmentation 
of 112 Sn, one can extract the relative chemical potentials between the two systems which appears to be robust 
with respect to sequential decay. Can these results be quantitatively understood within the framework of a 
statistical model? 

5. Might multifragmentation provide a superior environment for the production of exotic isotopes in certain regions 
of the N/Z-plane? Currently, abrasion/ablation models |Q-g| provide the preferred scenario for creating rare 
isotopes and experiments have focused on searching for rare fragments at projectile rapidities. The EPAX 
parameterization M, which is in common use for designing experiments, is built around such a picture. It is 
not clear to what degree multifragmentativc pictures can either complement or compete with abrasion/ablation 
models. 
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The canonical ensemble is employed in this study. The motivation for choosing the canonical ensemble, along 
with descriptions of the statistical fragmentation algorithm and the sequential decay calculation, appear in the next 
section. Isotope yie lds are primarily determined by the breakup temperature and sequential decay. These themes are 
introduced in Sec. |lll| . The importance of enforcing exact conservation of overall charge and baryon number with 
the canonical ensemble is discussed in Sec. Evl. The sensitivity to the size and neutron excess of the overall system is 
explored in Sec. [V| while Sec. VI contains a study of the sensitivities to several aspects of the model such as breakup 



density, level density and tunneling through the Coulomb barrier during the sequential decay process. Recently, Xu 
ct al. have measured yields of light fragments from the fragmentation of both 112 Sn and 124 Sn, and by taking ratios 
of isotope yi elds have determined the relative c hemic al potentials of the two systems [|| . This result is interpreted 



with in Sec. VII. A summary is provided in Sec. VIII 



II. MODELING STATISTICAL MULTIFRAGMENTATION AND THE SUBSEQUENT DECAY 

The initial statistical fragmentation of the system is calculated assuming the canonical ensemble. The choice of 



the canonical ensemble is explained in the following subsection. The recursive method described in Sec. [I B provides 



a pr obab ility for creating any state of any nuclide. The modeling of the subsequent sequential decay is described in 



Sec. [IC 



A. Choosing a statistical ensemble for fragmentation calculations 

Statistical calculations have proven remarkably successful in describing mass yields for multifragmentative processes. 
Such models assume a fixed volume where all distributions of N neutrons and Z protons into various fragments and 
excited states are considered. A rather wide variety of thermal models have been employed in the study of nuclear 
fragmentation, including grand canonical S, canonical lj and microcanonical |2| treatments. Additionally, isospin 
effects have been studied in lattice-gas models |j-[ll]] and in percolation |l2|;|l3[ . Since our goal is to study the yields 
of rarely produced isotopes which may be produced in future radioactive beam facilities as rarely as in one per 10 
events, we employ the methods recently promoted by Chase and Mekjian ]14|-[lq| which forego the need of using Monte 
Carlo methods. Majumder and Das Gupta have in fact studied Boron, Carbon and Nitrogen isotope production with a 
similar model to what is presented here |17| ] . The disadvantage of this method is that explicit interaction of fragments 
(beyond mean field or excluded volume effects) is outside the scope of the formalism. 

Furthermore, we choose to focus on the canonical ensemble rather than the microcanonical ensemble which would 
strictly constrain the energy of the sampled configurations. The canonical distribution considers all configurations of 
the same net N and net Z, weighted by the Boltzmann factor, e~" E . The mean energy and variance of the energies 
are given by derivatives of the canonical distribution, 

{H) = _d\og{n c ) 
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It has been shown that these systems undergo first-order phase transitions, which implies an infinite specific heat for 
infinite systems. However, for systems of size A w 200 the energy rises by approximately 4 MeV per nucleon as the 
temperature changes by one MeV in the range of T « 6 MeV. This results in a variance of the energy per nucleon, 

2T 

° E/A * vT (4) 



Thus canonical calculations effectively sample configurations within an excitation energy window of approximately 
1.0 MeV per nucleon, or less if the systems are much larger than A — 100, or if the temperature is away from the 
transition region. 

Unless an experiment gates on excitation energies to an accuracy better than 1.0 MeV per nucleon, there is little 
motivation to pursue a microcanonical treatment. Microcanonical treatments have been shown to give much different 
results than canonical treatments in the study of fluctuations and fragment multiplicity distributions [p.6[ , topics which 
will not be pursued in this study. 



Grand canonical treatments sum over configurations with different total Z and N. Since we are interested in the 
production of extremely neutron-rich or proton-rich fragments, it seems prudent to use the canonical ensemble. A 



comparison of grand canonical and canonical calculations is presented in section IV 



B. Fragmentation calculations with the canonical ensemble 

Our calculations are based upon the recursion relations of Chase and Mekjian which permit one to sum over all 
partitions of the system into different nuclear fragments. The recursion relation for the canonical distribution function 

is 

nz,jv(T) = X>-^« (5) 
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where the label c in the summation in the first line indicates one particular configuration, and Wi is the partition 
function of a specific nuclear species i, 
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internal levels j 

where m, is the mass and Zi, n, and a, are the charge neutron number and baryon number of the species i. The 
reduced volume, V Ie d — HotaJ — 2^4/pj,, roughly accounts for the overlap of the nuclei. Given the default breakup 
density of one sixth p$, the resulting reduced volume is four times the volume of a system at normal nuclear density. 
Calculating the partition functions of specific species u>i is straightforward given the levels and degeneracies of 
the nuclei. This was done for all fragments with a < 6. However, many of the nuclei which are generated in this 
approach have not yet been observed, and for almost all of them the level structure is not known. Given our lack of 
understanding of the ground state energies, let alone the excited state energies, we employ a liquid-drop treatment. 
We have chosen the finite-range liquid-drop model (FRLDM) |lq | as a means for generating ground state energies and 
have ignored the microscopic terms in the model which account for shell structure. The spectrum of excited states was 
generated by assuming a uniformly spaced assortment of single-particle states with spacing AE — a/^/A, where a is 
the level density parameter, chosen to be 10 MeV. The degeneracy of a state with excitation energy nAE was found 
by counting all ways to arrange particles and holes such that they summed to the desired excitation energy. For heavy 
nuclei, the separation AE becomes small, and when the separation fell below one MeV, the spectrum is interpolated 
onto a mesh of 1.0 MeV resolution. Finally, all mass- formula energies were modified to account for screening of the 
Coulomb potential, 

Seoul. -* Seoul. (1 - p/po), (7) 

where p/ po is the ratio of the density to nuclear matter density. 

Once the functions u>i(T) were calculated, the partition function £In.z(T) was then generated through the recursion 
relation, Eq. (pft, in less than a second of computer time. Subsequently, the populations of individual levels j of a 
species i were found, 
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where Ei t j and dij are the energy and degeneracy of the level. 



C. Modeling sequential decay 

The vast majority of the nuclear levels considered in these calculations are particle-unstable, including the ground 
states of those species outside the proton and neutron drip lines. After the initial yields were calculated, the subsequent 
decay was modeled by apportioning the weight of an unstable level into all the levels into which the nucleus might 
decay. Eight decay modes were considered: proton, neutron, deuteron, dineutron, diproton, t, 3 He and a. The decay 



weights were chosen according to Weisskopf arguments. Considering the decay from a state i to a state / via the 
emission of particle of type k with energy E). — Ef — Ei , 



Hi->f(E k ) ex mE k <Jkj-n, 



(9) 



where a^j^i is the capture cross section which we choose as the geometric cross section, n(Rf + Rk) 2 , scaled by the 
Coulomb correction, 



<7k,f->i = 7T(-R/ + Rk) 
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Here, E cou \ is the Coulomb barrier, 
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To account for particles tunneling through the Coulomb barrier, the s-wave contribution to tunneling through the 
barrier was estimated through the WKB approximation as 
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If the energy exceeded the Coulomb barrier, the capture cross section was taken as <j s = ir/k 2 . For all decays, both 
the classical form for the cross section in Eq. (jTo) and the s wavepiece described above were compared and the larger 
of the two cross sections was chosen for the decay rates in Eq. (0). 

Decays were calculated for all levels in all nuclei, beginning with the heaviest nuclei. For the decay of each level, 
the decay rate was calculated into every possible level energetically accessible through the eight decay modes listed 
previously. The weight associated with the decaying nucleus was then apportioned into all the states according to 
the rate for the decay into each state. The weights were also simultaneously added into the ground states of the 
eight nuclei representing the eight decay modes. Thus, the decaying process exactly preserved the initial N and Z of 
the original system. Given that heavier nuclei might have up to a thousand levels to decay with each level decaying 
into any of several hundreds of levels in a daughter nucleus, and given the existence of thousands of isotopes, the 
calculation of sequential decays often required a few hours. 



III. SENSITIVITY TO TEMPERATURE AND THE ROLE OF SEQUENTIAL DECAY 



Fragmentation is clearly sensitive to temperature, not only in the mass yields but also in the isotopic yields. Higher 
temperatures allow the system to produce nuclei with unfavorable symmetry energies more efficiently. However, 
higher temperatures also result in hotter nuclei that experience longer decay chains which in turn push the resultant 
residues closer to the valley of stability. The interplay of temperature and symmetry energy is discussed in the 
following subsection by studying the liquid drop model. Sec. Ill B presents results of the canonical calculation with 
and without sequential decay. By modeling the fragmentation and decay of a system with A = 200 and Z = 80, it 
is found that for lighter masses, a w 40, isotopes with highest n/a are best produced at temperatures near 5 MeV. 
Heavier neutron-rich fragments, a ss 100, are best produced at lower temperatures near 3 MeV, but such production 
is certainly in the domain of fission, and purely statistical calculations should not be taken too seriously. 



A. Insights from the liquid-drop model 



For a nuclear system at high excitation, where the initial nucleus is nearly completely dissolved into the constituent 
protons and neutrons, the neutron-to-proton ratio and the triton-to- 3 He ratios are roughly equal to the neutron-to- 
proton ratio of the decaying nucleus. At lower excitation, where most of the mass resides in large nuclear fragments, 



this is no longer true. In heavy nuclei, the relative populations for nuclides of fixed a are determined by the symmetry 
energy, temperature and Coulomb energy. Applying the grand canonical ensemble, one can make a liquid-drop 
estimate of the relative populations given the chemical potential, ^i = (/i„ — /i p )/2, related to the difference in 
neutron and proton number, 
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where £ is the symmetry term in the liquid-drop model. The Coulomb term also contributes to the symmetry energy. 
The parameters chosen here are b c = 0.7 MeV and £ = 23.4 MeV. A more complicated parameter set is employed in 
the FRLDM. However, for illustrative purposes, we consider the Bethe-Weizsacker liquid-drop model here. 

Ignoring the Coulomb term and completing the square in the exponential above allows calculation of the mean 
neutron excess for fragments of mass a. 



(n — z) = a- 
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The neutron to proton ratio is then, 
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If most nucleons reside in large fragments, the ratio (n — z) / a of a single species equals the ratio for the entire system. 
By considering the limit where (n — z) / a is small, one can derive the isospin amplification factor, 
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Thus, the isospin amplification factor is 2£/T, which exceeds 10 for temperatures less than 5 MeV. This ratio is the 
same for all mirror nuclei related by changing one proton to one neutron, e.g. t— 3 He and 15 0- 15 N. Large ratios of 
mirror nuclei have been observed in intermediate-energy heavy ion collisions 0]. It should bc noted that Serot and 
Miiller postulated similar ideas regarding isospin fractionation in terms of phase separation [pL9|| and that Samaddar 
and Das Gupta have made similar conclusions using the lattice gas model ]10| . 

The distribution described in Eq. (|l^) can be described by an offset and a Gaussian width. The offset is largely 
determined by the neutron fraction N/A of the composite system and does not depend strongly on the temperature 
unless the temperature becomes so high that protons, neutrons, tritons and 3 He fragments can absorb much of any 
initial excess of the neutron number. The width depends on the ratio of symmetry energy to the temperature. 
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Thus, although the ratio of yields of mirror nuclei tend to be largest for lower temperatures, higher temperatures 
result in the broadest widths and the greatest yields for nuclei with extreme neutron or proton imbalances. However, 
this conclusion will be significantly modified by sequential decay. 



B. Results for the fragmentation of A — 200, Z — 80 nuclei 

In this section we present results from the calculations described in Section O which combine the canonical ensemble 
with sequential decay. We will see that the expectations from the previous section, that higher temperatures produce 
more isotopes far away from the valley of stability, are qualitatively modified by the effects of sequential decay. 
Sequential decay will more strongly narrow the distribution when the temperature is higher and more decays occur, 
allowing the final residue to approach the valley of stability more closely. Combined with the production probability 
of producing a fragment of size a with any n — z, the optimum temperature for producing fragments with a large 
neutron excess is not trivially understood. 

Figure [l] displays the yield of a = 40 fragments as a function of n — z for the fragmentation of an A = 200, Z = 80 
system at several temperatures. The upper panel displays yields for the case where sequential decay is neglected. 



Clearly, higher temperatures create the best conditions for creating highly neutron-rich nuclei when sequential decay 
is neglected. The T = 4.5 MeV results are well described by the Gaussian resulting from the liquid-drop expression, 
Eq. (|13|). The magnitude and center of the Gaussian were adjusted to fit the peak but the width was determined by 
the liquid-drop parameters. 

Post-decay yields are displayed in the lower panel of Fig. |l|. Sequential decay strongly narrows the distribution, 
especially for higher temperatures. Figure g displays the yields of sulphur isotopes for the same system. From Figs. |l| 
and |2[ one concludes that the probability of producing rare neutron-rich isotopes is highest for temperatures near 5 
MeV for mass a « 40 fragments. 

Yields of a — 100 fragments are shown in Fig. |3| Since heavier fragments evaporate more particles, they are more 
affected by sequential decay. For heavier fragments, especially at temperatures less than 3 MeV, fission should provide 
the dominant production mechanism. Although the statistical description presented here does consider fission-like 
partitions of the system, both dynamics and shell structure play a significant role in fragmentation at low excitation 
where barriers play a pivotal role. For that reason, we confine our future discussion to the production of lighter 
fragments, a < A/3. 

IV. GRAND CANONICAL VS. CANONICAL ENSEMBLES 

The canonical ensembles employed in this study enforce exact conservation of neutron and proton number. In 
order to understand the importance of this constraint, we present corresponding calculations in the grand canonical 
ensemble. In the grand canonical ensemble densities are determined by chemical potentials, and the presence of a 
fragment of type i does not affect the possibility of observing a second fragment. A grand canonical ensemble supposes 
either an infinite system or one chemically connected to an infinite bath of both heat and particles. 

Performing calculations in the grand canonical ensemble begins, as in the canonical ensemble, by calculating par- 
tition functions as defined in in Eq. (|6j) . The yield of fragments of type i is then 

Y i =w i {T)e^ ni/T+ ^ ZiIT . (19) 

One must then find the chemical potentials fx p and /x„ that give the desired net neutron and proton numbers. In 
practice, this is accomplished numerically by applying Newton's method. For the results presented here, only nuclei 
of mass and charge less than that of the entire system were included in the sum. Neglecting higher masses only affects 
results at low temperatures. 

Yields of a = 40 fragments from grand canonical and canonical calculations are presented in Fig. for the frag- 
mentation of an A — 200, Z = 80 system, and a system of half that size, and are compared to results using the 
grand canonical ensemble. The yields are similar for the two ensembles for the fragmentation of the heavier system, 
but are significantly different for the lighter example. The reduced production of extremely neutron-rich nuclei from 
the fragmentation of the A = 100 system is easily understood as resulting from the absolute conservation of neutron 
number. In summary, use of the grand canonical ensemble is invalid for smaller systems, and even in the case of large 
systems, can result in overpredicting yields by an order of magnitude near the neutron drip lines. 

V. ISOTOPE PRODUCTION AS A FUNCTION OF SYSTEM SIZE AND CHARGE 

By fragmenting heavier systems, one is able to create hot systems with higher neutron fractions. By choosing 
larger systems, yields are less affected by the constraints of exactly conserving the net number of neutrons and 
protons. Although the qualitative trend is clear - larger and more neutron-rich systems are better suited to producing 
extremely neutron-rich fragments - quantitative dependencies might be rather surprising given the amplifications 
mentioned earlier. 

Yields for sulphur fragments are shown in Fig. for three systems: (A — 200, Z = 80), (A = 200, Z = 85), and 
(A = 100, Z = 40). Each system is chosen to fragment at a temperature of 4.5 MeV. The importance of choosing the 
most neutron-rich isotope is clearly apparent in the comparison of the Z = 80 and the Z = 85 systems. Here, the 
neutron fraction changes from 0.6 to 0.575. Furthermore, comparison with the smaller system at the same neutron 
fraction N/A = 0.6 demonstrates the modest enhancement of neutron-rich fragments involved in using larger systems. 

Given the sensitivity of the yields to changing the neutron and proton number of the fragmenting system, evaporation 
of neutrons during the early stages of the reaction should significantly reduce the production of rare neutron-rich 
isotopes. During the first 100 fm/c, Weisskopf arguments suggest that 20 to 30 neutrons, but only 5 to 10 protons, 
might evaporate from a Uranium nucleus heated to a temperature of 10 MeV. The neutron fraction of the residue might 
thus change from 0.614 to ~0.59 during the evaporation which significantly reduces the possibility to emit fragments 



near the neutron drip line. Fission yields are also affected by the number of evaporated pre-fission neutrons, which 
can be sizable, especially when the fission time scale is large. 

If one wishes to find reactions with the highest abundances of rare neutron-rich fragments, one should try to 
minimize the number of neutrons emitted in the pre-fragmentation stage of the collision. Reducing the number of 
neutrons emitted during the pre-fragmentation stage might be accomplished by choosing heavy-ion collisions as a 
means of producing multifragmcntcd systems. Since heavy ion collisions, especially mid-central collisions, retain some 
of the initial collective longitudinal motion of the beam, the expansion does not require as much time to accelerate 
from rest as it would in a pA collision. Furthermore, by colliding two very heavy nuclei, the global system would 
be relatively proton rich and not emit neutrons as preferentially, relative to protons, during the early stages of the 
collision. In fact, the mid-rapidity region, often referred to as the neck, might be even more neutron rich than the 
target or projectile [|20f . Of course, the disadvantage to using heavy- ion collisions is that lower beam energies are 
used and the produced fragments are not as forward focused which might make them more difficult to detect. If one 
wishes to produce fragments with charges greater than Calcium, z = 20, and energies less than 30a MeV, electron 
recombination and the resulting impure charge states makes detection and refocusing especially difficult. 

VI. SENSITIVITY TO BREAKUP DENSITY, LEVEL DENSITY AND SUB-BARRIER EVAPORATION 

Although the calculated yields of rare fragments depend on details of the modeling procedure, it is not obvious 
to what degree yields are sensitive to particular aspects of the modeling. One obvious parameter is the symmetry 
energy parameter. As the symmetry energy is uncertain by one or two MeV outside the region of measured nuclei the 
associated Boltzmann factors might change by a factor of two. Clearly, a lower symmetry energy will lead to broader 
isotope distributions. The sensitivity to other aspects of the modeling is not so trivial. In this section we consider 
three sources of uncertainty: the breakup density, the level density parameter and sub-barrier evaporation. 

Adjusting the breakup density affects the mass distributions as larger breakup volumes lead to preferential emission 
of many small fragments. To determine the sensitivity of the isospin distributions to the breakup density, the 
fragmentation of an A = 200, Z = 80 system was considered at three breakup densities, po/8, the default breakup 
density po/6, and po/4. The yield of fragments of mass a — 40 and of sulphur isotopes are displayed in Fig. @ assuming 
a breakup temperature of 4.5 MeV. The distributions of isotopes of fixed mass are fairly insensitive to the choice of 



breakup density. This is not surprising given the liquid-drop arguments presented in Sec. Ill A. The distribution 
of sulphur isotopes is broader for the higher breakup densities as would be expected since higher densities lead to 
increased production of heavier fragments. 

The level density parameter affects results with regards to the initial population of isotopes, the relative populations 
of levels of a given isotope, and in the dynamics of the sequential decay. The breakup of an A — 200, Z = 80 system 
is again considered, and the yields of mass 40 fragments and sulphur isotopes are displayed in Fig. [7| for three level 
density parameters, 8 MeV, the default value of 10 MeV, and 12 MeV. The distributions are noticeably broader for the 
higher level density parameters. Very little effect was evident in the pre-decay yields and thus this sensitivity derives 
from the dynamics of the sequential decay. This suggests the importance of improving the modeling the density of 
excited states as the level density is known to have a strong dependence on shell structure. 

The final aspect of the modeling considered here is barrier penetration, which plays an important role in sequential 
decay. Sequential decay, as seen in Figs. [I] and g, strongly narrows the isotope distributions as evaporation pushes 
the fragments towards the valley of stability. The Coulomb barrier also leads to preferential emission of neutrons 
relative to protons. In our calculations a fixed Coulomb barrier was imposed with an additional consideration of s- 
wave tunneling. To understand the importance of tunneling, calculations were performed with and without tunneling. 
Calculations were also performed with the Coulomb barrier completely absent. Results considering the fragmentation 
of the A = 200, Z = 80 system at a temperature of 4.5 MeV are displayed in Fig. g. Ignoring the tunneling increased 
the yields of proton-rich fragments but left the yield of neutron-rich fragments unchanged. Ignoring the Coulomb 
barrier altogether significantly reduced the population of proton-rich fragments and only slightly increased the yield 
of neutron-rich fragments. 

Interpretation of Fig. @ suggests that the yield of neutron-rich fragments is surprisingly insensitive to the imposition 
of the Coulomb barrier as the preference of emitting neutrons to protons is more strongly derived from the greater 
phase space available when decaying towards the valley of stability. If one were interested in the production of 
proton-rich fragments, there would be a strong motivation to improve the treatment of the Coulomb barrier. 



VII. ISOSPIN AMPLIFICATION AND BARYONIC CHEMICAL POTENTIALS 



The phenomenon of isospin amplification was presented in Sec. Ill A in the context of the liquid drop model and 
was also discussed in Sec. M. It was seen that relatively small neutron excesses could result in large increases in the 
yields of neutron- rich isotopes. Experimentally, it has been observed that yields are indeed quite sensitive to the 
neutron fraction of the fragmenting system. By creating a ratio of yields for different isotopes from the fragmentation 
of 112 Sn and dividing by the same ratio for the fragmentation of 124 Sn, one is able to extract the difference of the 
chemical potential for the two fragmenting systems. 

To better explain the ratio, one considers the population of a level i of an isotope of mass a in the grand canonical 
ensemble, 
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From the results of the previous sections, it is clear that the yield is significantly affected by sequential decay. However, 
studies using the MMMC model pllp| have revealed a surprising insensitivity of the following ratio to sequential decay, 
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where the 124 and 112 refer to the different isotopes of the overall system. This ratio should yield e ( n - z ) A t*™/ T ; where 
A/i„ is the difference of the neutron chemical potentials of the A = 124 and A = 112 systems. By creating a similar 
ratio using isotopes of the same neutron number, one is able to extract Afi p . 

To illustrate the insensitivity of the ratio to sequential decay, the logarithm of the ratio is shown both before and 
after sequential decay in Fig. |] for the fragmenting 112 Sn and 124 Sn systems at a temperature of 4.5 MeV. The slight 
non-linearity of the pre-decay results derives from the finite size of the system, i.e. a grand canonical calculation 
would yield a straight line. 

The resulting ratio also clearly depends on temperature and breakup density. To illustrate the temperature de- 
pendence calculations were performed for three temperatures, 3.0 MeV, 4.5 MeV and 6.0 MeV. The log of the ratio 
was in this case multiplied by the temperature so that one could focus on A// rather than A/i/T. The slopes weakly 
depend on the temperature as shown in the lower panel of Fig. HOL The sensitivity to the breakup density was also 
investigated and found to be rather weak. 

The upper panel of Fig. Wu displays a line derived from analysis of several experimentally measured isotope ratios 
J3|, where the ratios were observed to be surprisingly robust with respect to the choice of element and isotopes. The 
comparison with data is difficult to interpret. At the time in a collision where multifragmentation occurs, one would 
expect some of the extra neutrons in the 124 Sn system to have boiled off. Thus, it is not surprising that the calculation 
over-predicts A/i n as the calculation assumes the systems differ by 12 neutrons. The experimental result was also of 
a different effective size due to the fact that the mechanism for excitation was semi-peripheral heavy ion collisions. 
As the experiment involved 124 Sn + 124 Sn, 112 Sn + 124 Sn, and 112 Sn + 112 Sn collisions, the effective system size 
might have been either more or less than the size of a single Sn ion depending on the centrality of the collision and 
experimental acceptance. Nonetheless, several positive conclusions should be drawn from this analysis. First, the 
ratios are indeed robustly independent of sequential decay, even when one considers isotopes far along the isotope 
chain from stability. Secondly, the experimental results would seem to be consistent with model predictions if one 
were to assume that approximately a third of the extra neutrons gained from going from 112 Sn to 124 Sn evaporated 
prior to multifragmentation. 

VIII. SUMMARY AND CONCLUSIONS 

Statistical multifragmentation has been investigated in regards to the production of rare nuclides, particularly 
neutron-rich isotopes. Statistical calculations have been performed in the canonical ensemble, taking into account 
exact conservation of N and Z. By comparing to calculations in the grand canonical ensemble, particle- number 
constraints were seen to be non- negligible, especially for the fragmentation of smaller systems. Assuming a simple 
liquid-drop form for the ground-state energies and assuming a simplified form for the distribution and degeneracies of 
excited states, initial yields were calculated for every level of every possible isotope. Sequential evaporation has been 
modeled by calculating the ways to partition the initial weight of every particle-unstable level into any state to which 
it might decay. Eight modes of decay were considered. Since neither the calculation of the initial distribution nor the 



modeling of the sequential decay involved Monte Carlo procedures, it was possible to calculate yields for extremely 
rare isotopes. 

By considering yields of fragments of a fixed a as a function of n — z, two lessons were learned by considering 
a simple liquid drop model. First, the width of the initial distribution was determined largely by the ratio of the 
temperature to the symmetry term in the liquid-drop model. Secondly, isospin amplification, or isospin fractionation, 
could also be simply understood in terms of the same two quantities. 

Sequential decay dramatically altered the yields, pushing the yields towards the valley of stability. The effects 
were strongest for larger fragments and for higher temperatures. As the initial distributions for fragments of fixed 
a as a function of n — z were broadest at high temperature, inclusion of sequential decay shows that 5 MeV is 
the best temperature for creating rare neutron-rich fragments in the a — 40 region. For heavy fragments with 
a = 100, the optimum temperature was near three MeV which is more in the domain of fission than the domain of 
multifragmentation. 

In order to understand the sensitivity of the calculations to various aspects of the modeling, the breakup density, 
level density, system size and system neutron fraction were systematically varied. Results were weakly dependent on 
the breakup density and moderately sensitive to the level density. Larger systems lead to somewhat broader yields 
due to finite-size constraints. Yields were especially sensitive to small changes in the isospin composition of the 
overall system. Thus, choosing projectiles and targets with large neutron fractions, e.g. N/A for Uranium is 0.614 
and N/A for 124 Sn is 0.597, strongly increases the chances of creating isotopes near the neutron drip line. Finally, the 
sensitivity of the yields with respect to details of the evaporation was considered by adjusting the Coulomb barrier 
that leads to preferential emission of neutrons as opposed to protons. The yields of neutron-rich fragments were 
surprisingly insensitive to the barrier and changed almost imperceptibly when the Coulomb barrier was removed 
altogether. Proton-rich fragments were however quite sensitive to the details, as the tunneling allowed proton-rich 
fragments to return to the valley of stability. 

The ability of statistical models to explain isotope production is certainly of scientific interest in its own right. 
Additionally, one might also consider whether multifragmentation could offer a competitive means for creating rare 
isotopes. However, separating and focusing particles could be problematic. If one fragmented heavy nuclei, which have 
the broadest yields, the Coulomb forces between the fragment of interest and various parts of the residual system 
would spread the emission over a large kinematic region making separation or focusing of the produced particles 
difficult. This would be especially true for fragments produced at mid-rapidity. Thus, it is difficult to discern whether 
there is a pragmatic side to the isospin degree of freedom in multifragmentation. To date, experiments that measure 
isotope yields have either been designed to focus on projectile rapidities and thus ignore multifragmentative events, 
or have measured production of light elements, a < 20. If statistical descriptions are shown to accurately describe 
isotope production for a broader range of nuclides, it might warrant serious consideration of designing an apparatus 
to capture and identify very rare isotopes in a very different environment than projectile fragmentation. 
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FIG. 1. Yields of mass 40 fragments are shown for the fragmentation of an A — 200, Z = 80 system at a variety of 
temperatures. Yields are shown both for the case where sequential decay is included (lower panel) and neglected (upper 
panel). The pre-decay yields at 4.5 MeV are well described by a Gaussian (line) where the width is determined by liquid drop 
parameters described in Eq. (|l3|). Pre-decay yield curves are broadest for higher temperatures, but post-decay yields of very 
neutron rich fragments are largest for temperatures near 5 MeV. 
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FIG. 2. Sulphur yields from the fragmentation of an A — 200, Z = 80 system are displayed for a variety of temperatures. 
Yields are shown both for the case where sequential decay is included (lower panel) and neglected (upper panel). As in Fig. [fj 
the best chance of producing isotopes near the neutron drop line occurs for temperatures near 5 MeV. 
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FIG. 3. The same as Fig. hi except for mass 100 fragments. Post-decay yields for extremely neutron-rich fragments in this 
mass range are largest for temperatures near 3 MeV. 
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FIG. 4. The importance of enforcing baryon number and charge conservation in the statistical treatment is illustrated 
by comparing yields of sulphur isotopes for canonical and grand canonical calculations. As expected, fragmentation of the 
A = 200, Z — 80 system is better described by the grand canonical calculation than fragmentation of the lighter system, 
A = 100, Z = 40. 
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FIG. 5. Yields of mass 40 fragments (upper panel) and sulphur isotopes (lower panel) are shown for a variety of systems 
fragmenting at a temperature of 4.5 MeV. Yields of neutron-rich isotopes were moderately decreased when considering a smaller 
overall system of the same neutron fraction. For the A = 200 systems yields were substantially reduced by lowering the net 
number of neutrons from 120 to 115. 
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FIG. 6. Fragmentation of an A — 200, Z = 80 system at a temperature of 4.5 MeV is calculated for three different breakup 
densities. Yields are shown for mass 40 fragments as a function of n — z (upper panel) and for sulphur fragments as a function of 
mass (lower panel). The breakup density affects the mass distribution more strongly than the relative probabilities for isotopes 
of fixed mass. 
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FIG. 7. The yields of mass 40 isotopes (upper panel) and sulphur isotopes (lower panel) are shown to be sensitive to the 
choice of level density parameter a. Higher level density parameters lead to broader isotope distributions due to the dynamics 
of sequential decay. Again, calculations concerned an A = 200, Z — 80 system at a temperature of 4.5 MeV. 



I 1 

o 
^- 

32 

CD 

*>- 

i i 

o 

cn 
O 







-5 



-10 - 



-15 




-10 







10 



20 



n— z 



D 

CD 

• — 

>- 



CD 
O 




I ' I l_ 

a No Barrier - 

o No tunnelingj 

n Default I 



1 



30 40 50 60 



a 



70 



FIG. 8. The sensitivity of isotope yields to details of evaporation criteria are studied by considering the yields of sulphur 
isotopes from the fragmentation of an A = 200, Z = 80 system at a temperature of 4.5 MeV. Results from three cases are 
displayed. The evaporative description in the default calculation (squares) assumes a Coulomb barrier with an s-wave tunneling 
contribution described in Eq. ([13). Calculations were repeated with no tunneling (diamonds) and with the barrier eliminated 
(triangles) , making the emission of neutrons and protons equal. Yields near the neutron drip line were insensitive to details of 
barrier, but yields near the proton drip line were indeed sensitive to the barrier and to the tunneling. 
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FIG. 9. Ratios of Carbon isotope yields for the fragmentation of 12 Sn are divided by the same ratios for the fragmentation 
of 112 Sn to obtain a quantity which should reflect exp{(n — z)A/i n /T}, where A/i„ is the increase of the neutron chemical 
potential resulting from using a heavier tin isotope. The logarithm of that ratio, which is presented here, rises linearly with a 
as expected. The pre-decay results (squares) and the post-decay results (squares) are remarkably similar. 
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FIG. 10. In the top panel the same quantity as in Fig. M is displayed for calculations at three temperatures, 3.0, 4.5 
and 6.0 MeV. The solid line represents the experimentally extracted value of A/j,„ from reference M. The consistency of the 
experimental results with the T = 3 MeV results may be misleading as the calculations ignored the fact that some of the extra 
neutrons associated with 124 Sn could have evaporated before the onset of multifragmentation. The same results are displayed 
in the lower panel, only multiplied by the temperature. This illustrates that Afi n is fairly insensitive to the temperature. 



